Nonuniversality of weak synchronization in chaotic systems 
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We show that the separate properties of weak synchroniza- 
tion (WS) and strong synchronization (SS), reported recently 
by Pyragas [K. Pyragas, Phys. Rev. E, 54, R4508 (1996)], 
in unidirectionally coupled chaotic systems, are not generally 
distinct properties of such systems. In particular, we find 
analytically for the tent map and numerically for some pa- 
rameters of the circle map that the transition to WS and SS 
coincide. 



Chaotic systems, by definition, are characterized by 
strong sensitivity to the initial conditions. Thus, in a gen- 
eral situation, one cannot synchronize two chaotic sys- 
tems, since in a practical situation it is impossible to start 
the evolution of the two systems with exactly the same 
initial conditions. However, one can synchronize subsys- 
tems of a chaotic system that have a chaotic output, but 
are intrinsically stable. This was shown by Pecora and 
Carroll |l[] . More specifically, Pecora and Carroll studied 
dynamical systems of the type u = g(u, w), w = h(u, w) 
and showed that variable w' governed by ib' = h(u, w') 
can synchronize with w if the sub-Liapunov exponents of 
the driven subsystem w' are all negative (implying that 
the subsystem w' is stable). The sub-Liapunov expo- 
nents they defined depend on the Jacobian matrix of the 
w subsystem, taking derivatives with respect to w only. 

In a recent letter, Pyragas [B studied one-dimensional 
chaotic systems, governed by a function /, coupled uni- 
directionally in the following way 



X n +1 = f(x„) 

Vn+1 = f(Vn) ~ k[f(Vn) ~ f(Xn)] 
Zn+l = f(z n ) - k[f( z n) ~ f(x n )]. 



(1) 



He found numerically that in a specific system, namely, 
the logistic map, synchronization between the variables 
x, y and z, as k is increased, occurs in two stages. In 
the first stage, y synchronizes with z and not with x, 
starting at a critical value k w . He called this weak syn- 
chronization (WS). The second stage of synchronization 
starts at k = k s , with k s > k Wy and is characterized by 
the synchronization of x, y and z, which he called strong 
synchronization (SS). Note that in above equation the 
case k = 1 represents synchronization in a trivial way, 
since this leads to y n +i = 2„+i = x n +i- So, for this kind 
of coupling k is taken in the interval [0, 1). 

The conditions under which WS and SS occur were 
determined by two Liapunov exponents, that is, the con- 
ditional Liapunov exponent, 



1 - 

= ln(l-fc) + iv lim ) -^l n |/'(y n )|, 



(2) 



defining the stability of the invariant manifold y = z, 
and the transverse Liapunov exponent of the invariant 
manifold x = y, 



1 N 

A =ln(l-fc)+ lim -£ln|/'(x n )|. 

N—>oo iv ,£ — ' 

n=l 



(3) 



In the logistic map, Pyragas found that X R becomes zero 
at two characteristic values of the coupling strength, k w 
and k s , corresponding to the threshold of WS and SS, 
respectively. In the region of weak synchronization (k w < 
k < k s ), \ R < and Ao > 0. Strong synchronization 
occurs for k > k s , where these two Liapunov exponents 
coincide, i.e., Ao = X R - 

We studied the phenomenon of chaotic synchroniza- 
tion in other one-dimensional maps, coupled also accord- 
ing to Eq. (0). We found that the phenomenon of weak 
synchronization is not always found in such systems. 

We start by showing what is the relationship between 
the exponents X R and Ao, defined by Pyragas, and the Li- 
apunov exponents of the three dimensional system (with 
dynamical variables x, y and z) governed by Eq. (jl]). We 
will call this the global system, which has the following 
Jacobian matrix at a single position along the orbit: 



f'(Xn) 
Jn = ( kf'(x n ) [1 - k]f(y n ) 

kf(x n ) 






[1 - k]f'(z n ) 



(4) 



The Liapunov exponent of this system is found by calcu- 
lating the eigenvalues of the matrix that consists of the 
product of the Jacobian matrices along a given orbit. It 
turns out that, because of the symmetry of this particu- 
lar matrix, the product J of the Jacobian matrices J n is 
of the type 



J 



an 
a 2 i a 2 2 
a3i a 33 



(5) 



where an, 021, a 2 2, 031 and 033 are in principle different 
from zero. The eigenvalues of this matrix are 



N 



Ai = an 



n /'^)' 



n=l 



N 



A 2 = aaa = IJ[l-*]/ , (lfrO, 



(6) 



(7) 



1 



A 3 = a 33 = l[[l-k]f'(z n ). 



(8) 



Consequently, in this way of coupling, the eigenvalues of 
the system are easily found, and each one is a function of 
a single variable. The Liapunov exponents are therefore 



1 N 

n=l 



A 2 = ln(l - k) 



1 



N 



71=1 



A 3 = ln(l -k)+ Urn 



N- 
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n=l 



ln|/'(^)|- 



(9) 



(10) 



(11) 



If the initial values of y and z are in the same basin of 
attraction, then A3 = A2, because the parameters of the 
maps are the same. The Liapunov exponent X R is equal 
to one of the Liapunov exponents of the global system, 
namely, it is equal to A2. Thus, this Liapunov exponent 
has a clear physical importance even when y and z are not 
synchronized. On the other hand, Ao is not a Liapunov 
exponent of the global system in the region where x and 
y are not synchronized (weak synchronization). However, 
Ao and Ai are related to each other via a trivial additive 
term, namely, A = ln(l — k) + Ai. It is obvious that in 
the region of SS, A = X R = A 2 . 

We found in fl that the loss of synchronization be- 
tween subsystems of a dynamical system coupled accord- 
ing to Pecora and Carroll's coupling is characterized by 
the the Liapunov exponents of the global system. That is, 
when one of such Liapunov exponents crosses zero (in the 
positive direction) the two identical subsystems will lose 
synchrony. When this occurs, a transition from chaos to 
hyperchaos Q takes place. Here we find a similar phe- 
nomenon. The transition from SS to WS corresponds 
to the transition from chaos to hyperchaos in the global 
system. 

Our next step is to show that WS does not necessarily 
precede SS. We show this analytically in the tent map, 
and later numerically in the circle map. 

The tent map is defined as 



/(*) = 



ax, 
a(l — x) 



if < x < 1/2, 
if 1/2 < x < 1. 



(12) 



The tent map has a period-one orbit when < a < 1, 
with fixed point x* = 0. For 1 < a < 2 the map has a 
chaotic orbit. For a > 2 the orbit diverges. (Note that, to 
avoid divergencies when < a < 2, the initial conditions 
must be in the interval [0, 1]). Since \f'(x n )\ = a in any 
point of the orbit, we have, using Eqs. (H)-(O) 



Ai = In a, 



A 2 = A 3 = m(l - k) + In a. 



(13) 



(14) 



For a fixed a, we find that A 2 decreases monotonically 
as k is increased from k — 0. Consequently, there is no 
region of WS in this map. For this map, synchronization 
between x and y and z occurs simultaneously, and only 
what is called SS is seen. The transition where SS occurs 
is determined by A 2 < 0, which gives k s = 1 — 1/a. 

Next, we consider coupled circle maps according to 
Eq. (pi), where 



lu sin(27rx) 

27T 



(15) 



Here we find regions of the parameter space where WS is 
not observed. We show in Fig. |l](a) an example for this, 
where 6 = 6 and lu = 0.44. The solid line in that figure 
represents X R , which is equal to A 2 , and the dashed line 
represents Ao, which is equal to log(l — k) + Ai. We do 
not find a region of k where X R < and Ao > 0, and 
consequently no WS is seen. 

For the circle map (and also in the logistic map) we ob- 
served interesting phenomena when the driving variable 
x is in one of the periodic windows of the chaotic band. 
There, we see regions of k where y and z are chaotic 
(A 2 = X R > 0) even when x is periodic (Ai < and 
Ao < 0). This is shown in Fig. 0(b), where b — 4 and 
lu = 0.4. 

In Pyragas also studied the coupling of the 

Rossler's and Lorenz's systems according to 



±1 = —a[x 2 + X3], 

x 2 = a[xi + 0.2x2], 

X3 = a[0.2 + x 3 (xi - 5.7)], 

Vi = 10(-yi +2/ 2 ), 

y 2 = 28yi - y 2 - ym + kx 2 , 

2/3 = 2/12/2 — 8/3y 3 , 



(16) 



(17) 



where a — 6. Here, we also can understand the phe- 
nomenon of weak synchronization reported by Pyragas 
for this coupled system, by considering it as a single 
global system of six variables. We calculated all the Li- 
apunov exponents of the global system, and found that 
at k w 6.6 one of the Liapunov exponents change sign. 
Beyond this critical value of k, synchronization between 
two copies of the Lorenz system [Eq. (0)] with different 
initial conditions (in the same basin of attraction) will 
occur. That is, WS will be seen. No other change of sign 
was found in the Liapunov exponents as we increased k to 
k = 200. Next, we calculated the information dimension 
Di of the attractor of the global system, assuming that 
the Kaplan and Yorke [|| conjecture holds. Contrarily to 
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the numerical results reported in Eq| , we find numerically 
that the dimension of the global attractor does not con- 
verge to the dimension of the driving system [Eq. ([l6|)]. 
The convergence of these two dimensions was considered 
in 0] as the characterization of SS in this system. As 
Fig. EJ shows, for k > 50, Di for the global attractor 
(solid line) is approximately constant, but different from 
the information dimension of the driving system (dashed 
line). For the global system we find that Di = 2.16 when 
k > 50 and for the driving system Di — 2.01 (with the 
last digit being uncertain in both cases). Consequently, 
although we observed the property of WS in the cou- 
pled system governed by Eqs. ([If]) and (|l7|), we have not 
identified a regime of SS, as characterized in 

In summary, we have found that the property of weak 
synchronization in maps coupled according to Eq. (|l|) 
is a only a particular property of some systems. We 
have found that it does not hold in at least two sys- 
tems, namely, the circle map and the tent map. We have 
also clarified the relationship between X R and Ao and the 
Liapunov exponents of the global system. Finally, we 
have verified that strong synchronization as defined in 
[EJ is not observed in the coupled system governed by 
Eqs. (U) and |i]]). 



FIG. 2. Information dimension Di as a function of k for 
the global attractor [Eqs. (|l6|) and (0)] (solid line) and for 
the driving system [Eqs. (Jlq)] (dashed line). The error bars 
for Di are smaller than the diamond symbol. Since as k in- 
creases, the variables yi change faster than Xi we decrease the 
integration time step At according to At = 0.02/(fc + 1). The 
initial conditions used were X\ = 0.1, X2 = 0.2, x% — 0.3, 
2/i = 0.4, 2/2 = 0.5 and y 3 = 0.6. 



* Electronic address: mariav@eecs.berkeley.edu 
[1] L. M. Pecora and T. L. Carroll, Phys. Rev. Lett. 64, 821 
(1990). . 

[2] K. Pyragas, Phys. Rev. E 54, R4508 (1996). 

[3] M. de Sousa Vieira, A. J. Lichtenberg and M. A. Lieber- 
man, Phys. Rev. A 46, R7359 (1992). 

[4] O. E. Rossler, Notes Appl. Math. 17, 141 (1979); O. E. 
Rossler, Phys. Lett. 71A, 155 (1979). 

[5] J. Kaplan and J. Yorke, in Functional Differential Equa- 
tions and the Approximation of Fixed Points, edited by 
H.-O Peitgen and H.-O Walther, Springer Lecture Notes 
in Mathematics, vol. 730 (Springer- Verlag, Berlin, 1979), 
p. 204. 



FIG. 1. Liapunov exponents (in Pyragas' notation) X R 
(solid line) and Ao (dashed line) for coupled circle maps, with 
(a) b = 6 and lu = 0.44, and (b) b = 4 and lu = 0.4. We used 
N = 30, 000 in Eqs. (El, (|), and neglected a transient of 3,000 
iterations. The initial conditions were x = 0.1, y = 0.2 and 
z = 0.3. 
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Fig. 1(a) 
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